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Experimental observations of an almost symmetric electronic circuit show complicated sequences 
of bifurcations. These results are discussed in the light of a theory of imperfect global bifurcations. 
It is shown that much of the dynamics observed in the circuit can be understood by reference 
to imperfect homoclinic bifurcations without constructing an explicit mathematical model of the 
system. 
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I. INTRODUCTION 

The role of symmetries in determining the behaviour of 
nonlinear physical systems can be crucial. Reflection (or 
Z2) symmetry is relevant to a wide range of experiments, 
and in such a system a pair of stable solutions may be 
created by a supercritical pitchfork bifurcation as a pa- 
rameter is varied. These new states break the original 
symmetry, but are symmetric images of each other. Of 
course, perfect symmetry is never achievable in any phys- 
ical system so in practice the bifurcation may become 
disconnected having one branch which varies monotoni- 
cally with the parameter and a second which arises by a 
saddle-node bifurcation. This is most easily modelled by 
adding an imperfection term as a constant in the model 
normal form and this appears to work well in describ- 
ing the local bifurcation structure. However, a physical 
system will typically contain many sources for this imper- 
fection and some of them may be high-dimensional in na- 
ture. Therefore, it is reasonable to ask whether a model 
with a single imperfection term provides a good repre- 
sentation of the system far from the bifurcation point. 
Specifically, we are interested here in the effects of this 
local modelling on the global dynamics which result from 
homoclinic bifurcations. 

Our investigation is concerned with a class of global bi- 
furcations involving homoclinic orbits, i.e. orbits which 
tend to a stationary point of the model flow in both for- 
wards and backwards time. Typically, the existence of a 
homoclinic orbit is not a persistent property of a differen- 
tial equation, but they occur on lines in two-parameter 
families (technically, they are codimension one bifurca- 
tions) . In the absence of symmetry, the net effect of such 
bifurcations is to create or destroy a periodic orbit, whose 
period tend to infinity at the bifurcation point. This may 
happen in one of two ways: one-sided or two-sided. In the 
one-sided case, the orbit apporaches the bifurcation point 
from one side of the bifurcation point as its period tends 
to infinity. In the two-sided case, such as the Shil'nikov 
case , the locus of the orbit in parameter space oscil- 
lates about the bifurcation value creating the so-called 



'Shil'nikov wiggle' as the period of the orbit tends to in- 
finity. Moreover, there are period-doubling and reverse 
period-doubling bifurcations of the orbit together with 
more complicated homoclinic bifucations. This sequence 
of events has been reported previously |^ in an experi- 
mental and theoretical study of a modified van der Pol 
oscillator, and in a wide variety of other experiments in- 
cluding Taylor-Couette flows Q, optics, [El ^, chemi- 
cal oscillators |Q ^ and liquid crystal flows . 

In the presence of simple symmetries, homoclinic bi- 
furcations may involve two or more homoclinic orbits. In 
the simplest cases the net effect is to destroy a pair of pe- 
riodic orbits which are the image of each other under the 
symmetry and create a single symmetric branch of peri- 
odic orbits. These symmetric periodic orbits cannot un- 
dergo period-doubling bifurcations in the two sided case. 
The period-doubling and reverse period-doubling bifur- 
cations on branches of the symmetric orbit are replaced 
by an initial symmetry-breaking (or reverse symmetry- 
breaking) bifurcation. The asymmetric orbits created in 
this way may, of course, be involved in period-doubling 
bifurcations. This distinction will be useful in the inter- 
pretation of the bifurcations observed below. 

Whilst the effect of small symmetry-breaking terms on 
the bifurcations of stationary solutions has a long his- 
tory (the imperfection theory of Golubitsky and Schaef- 
fer |l^) there appears to have been no system- 

atic attempt to describe the equivalent modifications of 
global bifurcations (although see for a special 

case). Our aim here is to provide the foundations for 
such an approach. We reconsider the experimental elec- 
tronic oscillator |^ which exhibits a variety of almost 
symmetric global bifurcations and show how many fea- 
tures observed in the experiments may be explained by 
reinterpreting some results on codimension two homo- 
clinic bifurcations so as to obtain a general imperfection 
theory for homoclinic bifurcations. These results neces- 
sarily involve non-stationary solutions, and so are likely 
to be applicable and observable in many more interesting 
situtations. 

The experiments were carried out using a van der Pol 
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oscillator. The bifurcation structure of this system has 
been investigated in detail previously Q but with the 
implicit assumption of symmetry. It is the aim of the 
present study to investigate the global dynamics of the 
circuit and relate the observations to modern ideas on 
gluing bifurcations where the mathematical abstraction 
of perfect symmetry is relaxed. 

II. EXPERIMENT I 
A. The electronic oscillator 

The experimental study was performed using a van 
der Pol oscillator circuit, the details of which are given 
in Healey et al. It comprises an autonomous LCR 
oscillator with two nonlinear conductances in the feed- 
back circuit. Precise variation of the two parameters 
which control the behaviour of the system was provided 
by switchable decades resistance boxes. By this means 
determination of the bifurcation structure to a relative 
accuracy of better than 0.1% was possible. The two pa- 
rameters are denoted by ai,/?i and they are nondimen- 
sionalised forms of the resistances Ri , R2 which control 
the nonlinear elements. Details of the nondimensionali- 
sation are given in Healey et al. Q. 

The principle set of observations were made using 
an oscilloscope. Steady bifurcations were observed as 
changes in the level of the d.c. output. On the other 
hand dynamical states were best monitored as Lissajous 
figures formed from a combination of signals measured 
over the nonlinear elements. In this way, limit cycles, 
period doubling sequences, chaos etc. were readily dis- 
played. Time-series were also recorded and stored on a 
computer via a 12-bit A/D for further processing. This 
included phase portrait analysis using the method of de- 
lay coordinates. 

The inductor used in the present circuit is 1.5269H 
compared with 1.78H used by Healey et al. ||^. This 
causes a shift of the bifurcation points relative to those 
previously reported, though the bifurcation structure 
remains qualitatively the same. The imperfections in 
the circuit are tiny and the resulting disconnections are 
equally small. They arise from a variety of sources but 
we will refer to them throughout as a single imperfection. 

B. Bifurcation set 

The stability diagram for the electronic circuit is shown 
in Figure |l|. The overall structure shows lines of steady 
and dynamic bifurcations all meeting at the top right 
hand corner of the figure which is a codimension-2 point. 
The dynamic bifurcations (Hopf and homoclinic) are 
pairs of lines superposed and separated by the imperfec- 
tions in the circuit. This effect is very small and cannot 
be resolved on the scale of the figure but, as we will show 
below, it has a significant effect on the global dynamics. 
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FIG. 1: Experimental bifurcation set in the q:i,/3i plane. SN 
denotes the path of saddle-node bifurcations, 'Hopf the Hopf 
bifurcations to simple oscillations and 'Hom' the gluing bi- 
furcations. The paramater region denoted by 'P2' is where 
forward and reverse period doubling is observed on the asym- 
metric orbits. 

In the parameter range of interest, a perfectly sym- 
metric system would have a trivial zero volts fixed point 
which would lose stability to a pair of non-zero d.c. states 
at a supercritical pitchfork bifurcation. As expected, in 
the experiment we see that this bifurcation is discon- 
nected to form a continuously connected state and a sep- 
arate solution branch which is terminated at its lower end 
by a saddle-node bifurcation denoted by SN in Figure [|. 
The stable non-trivial asymmetric d.c. states both be- 
come time-dependent via Hopf bifurcations; one on each 
branch. The imperfection in the circuit is very small, 
so the loci of these bifurcations almost coincide and are 
marked 'Hopf in Figure |^. The two asymmetric limit 
cycles which arise at the Hopf bifurcations appear to 
glue together leading to a large symmetric periodic or- 
bit. This transition is denoted by the line marked 'Horn' 
in Figure |^ and will be discussed in detail below. This 
symmetric limit cycle undergoes different types of bifur- 
cation including symmetry-breaking and period doubling 
and may also become chaotic. Finally, within the oscil- 
latory regime forward and reverse period doubling se- 
quences have been observed and these can be related to 
the Shil'nikov wiggle as shown by Healey et al The 
boundaries of this region are denoted by P2 in Figure 



C. Imperfect gluing bifurcation 

We first examine the influence of the imperfection on 
the gluing bifurcation which occurs when the two asym- 
metric limit cycles join without the presence of compli- 
cated dynamics. We chose /3i sufficiently large (/3i > 0.59 
approximately) and ai close to /3i so that the chaos which 
arises from period-doubling sequences on a Shil'nikov 
wiggle is avoided and the dynamics is almost planar. We 
present a 'typical' set of results for the orbit structure 
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FIG. 2: Oscillation period of different periodic orbits at (i\ — 
0.6000 plotted as a function of ai. '1' and '0' denote the orbits 
on the asymmetric branches and '10', '01' are the glued orbits. 




form two homoclinic orbits, biasymptotic to the station- 
ary point. 

As is clear from Figure and as should be expected 
of a real physical system, the oscillator is not perfectly 
symmetric. Hence it is not surprising that the pair of ho- 
moclinic orbits which exist at a single parameter value in 
the symmetric system seem to occur at different param- 
eter values in the oscillator. The results shown in Figure 
^ also suggest that there is a third homoclinic bifurcation 
- the bifurcation which creates the large amplitude '10' 
periodic orbit. 

It can be seen in Figure || that the period of both the 
small asymmetric orbits '1' and '0' increases as ai in- 
creases and they finally lose stability and jump to the 
'10' orbit at ai « 0.6045 i.e. where the graphs of the 
variation of period are almost vertical. Moreover, the '0' 
orbit remains stable for slightly higher values of ai than 
the '1' orbit, emphasising that the two orbits are not the 
images of each other under the symmetry. It should be 
noted that the '1' orbit results from a Hopf bifurcation 
on the monotonic branch of the disconnected pitchfork 
bifurcation. Therefore it loses stability before the '0' or- 
bit. This is precisely what is predicted by the addition of 
a constant term to the normal form. The orbits shown 
in Figure || all coexist at ai ~ 0.6041 and are typical 
examples of the limit cycles involved in this gluing bifur- 
cation. The fact that they can all coexist explains why 
hysteresis can be observed in the experiments. 

There are three features in Figure ^ which we will seek 
to explain theoretically in the next section: the break up 
of the gluing bifurcation, hysteresis, and also the extra 
bifurcations evident at larger values of ai. Before de- 
scribing the theory we shall look at this latter sequence 
of bifurcations in more detail. 



FIG. 3: Phase portraits of coexisiting asymmetric (1,0) and 
symmetric (10) periodic orbits at ai — 0.6041 and /3i = 
0.6000 



of the oscillator in this regime in Figures g and |^ which 
were taken at /3i = 0.6000. Figure || shows the period 
of the various simple orbits observed as a function of the 
parameter ai, and Figure |^ shows the form of the cor- 
responding orbits - the two small asymmetric orbits are 
labelled by '1' and '0' respectively, and the large ampli- 
tude orbit is labelled by '10', for reasons which will be 
explained below. 

If the electronic oscillator were symmetric then the de- 
velopment of the orbits shown in Figure ^for ai = 0.6041 
would have a simple explanation in terms of gluing bifur- 
cations : two periodic orbits which are the symmetric 
image of each other approach a stationary point and are 
'glued together' to form the single symmetric orbit with 
code '10'. At the bifurcation the two smaller periodic or- 
bits touch at the stationary point, i.e. they are no longer 
periodic (their period has diverged to infinity) and they 



D. Symmetry- breaking bifurcation of large 
periodic orbit 

It is known that symmetric systems cannot undergo 
period-doubling sequences directly [l^ but must first 
break their symmetry. Hence, we would expect the large 
symmetric orbit formed by the gluing of the two asym- 
metric ones will suffer a symmetry breaking bifurcation, 
as predicted for the symmetric Shil'nikov wiggle [ pT[ . 
This was observed at /3i — 0.6000 for ai just above 
0.6059. The bifurcation was detected by measuring the 
mean voltage averaged over 150 periods of the oscillation 
and plotting this as a function of ai. The resulting bifur- 
cation diagram is shown in Figure ^ where we see that it 
has the form of a disconnected pitchfork. This diagram 
explains the creation of the orbit labelled '01' in Figure ||. 
Note that the original '10' orbit has a larger period but 
smaller < Vi > than the newly created '01' orbit and so 
the branches in Figures ^ andjj are apparently reversed. 
Two typical asymmetric orbits on respective branches are 
shown in Figure |for (ai, /3i) = (0.6067, 06000). It is ev- 
ident that the '10' orbit on the connected branch displays 
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FIG. 4; Bifurcation diagram of symmetry-breaking bifurca- 
tion of periodic orbits at Pi = 0.6000. The mean of Vi over 
5000 data points is plotted. 




FIG. 5: Phase space portrait of coexisiting "large" periodic 
orbits 10 and 01 at /3i = 0.6000 and qi = 0.6067 



strong variation in period for ai > 0.6065 and then insta- 
bility. However, the '01' orbit is virtually constant over 
this range. Both orbits then show reduction in period 
for high ai values. Each orbit undergoes period doubling 
sequences to chaos for ai values greater than the range 
displayed in Figure ^. The extra complications of period 
doubling and instability are topics for future research. 



III. THEORY 

It is natural to think of the bifurcations observed in 
the system in terms of two parameters. One of these, 
H say, is the parameter of the (fictional) symmetric sys- 
tem which has a gluing bifurcation as described in sec- 
tion lie. The second parameter, e say, is a measure of 
how far the oscillator is from being perfectly symmetric, 
i.e. it is some measure of imperfection with e = cor- 



responding to the perfectly symmetric system. Just as 
the standard imperfection theory for the bifurcations of 
stationary points ||l2| allows one to describe the effect of 
asymmetry in terms of /i and e, our aim here is to give an 
analogous description for general global bifurcations. We 
note that this is in the spirit of the work of Glendinning 
[Q and Cox jl3j for the particular case of Lorenz-like 
bifurcations. 



A. The basic picture 

Suppose that (/i, e) = (0, 0) denotes the point in pa- 
rameter space at which there are two symmetrically re- 
lated homoclinic orbits. Consider either one of these or- 
bits. Since the existence of homoclinic orbits is codi- 
mension one, there will be a curve in parameter space 
through (0, 0) on which systems have a homoclinic or- 
bit which is a continuation of the given orbit. Thus, for 
typical two-parameter families of systems, there will be 
two curves of homoclinic orbits in parameter space, Gq 
and Gi say, which intersect at the origin and which do 
not intersect the line e = again locally. The curve Go 
(respectively Gi) is the locus of a homoclinic bifurcation 
which creates or destroys the periodic orbit with code 
(respectively, 1). The one-parameter families of nearly 
symmetric systems such as the example considered in the 
previous section would then correspond to some curve in 
this two parameter space which has, for example, e > 
and which passes close to (/i, e) = (0, 0). Such a curve will 
intersect both Gq and Gi , but at different parameter val- 
ues, so there will be two simple homoclinic bifurcations 
at nearby parameter values on such a path. 

The intersection of the loci of two homoclinic bifur- 
cations (each to the same stationary point) is a codi- 
mension two phenomenon which has been studied by a 
number of authors |l|, |l|, |l|, |20[ |l|, ||, H, ||. The 
most important feature which all these bifurcations have 
in common is that at least two other curves of homoclinic 
orbits emanate from the intersection of Go and Gi , one in 
e > 0, labelled Gio, and the other in e < labelled Gqi- 
The labelling describes the order (in time) that the orbit 
passes through neighbourhoods of the basic homoclinic 
orbits. These homoclinic orbits are precisely the bifurca- 
tions needed to destroy or create (asymmetric) periodic 
orbits with code '10' or '01'. Thus a typical path close 
to e = will intersect Go, Gi and one of the curves Gqi 
or GiQ. This explains the third homoclinic bifurcation 
observed in Figure |^. Roughly speaking, the difference 
between orbits created by paths crossing Gio and those 
created by crossing Gqi is the difference between the or- 
bits shown in Figure ^. 

The details of the two-parameter bifurcation plane 
close to the intersection of Go and Gi depends upon the 
nature of the stationary point, the configuration of the 
homoclinic orbits and a measure of the amount of twist- 
ing of solutions about these orbits. The nature of the 
stationary point is determined by the eigenvalues of the 
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FIG. 6: The two parameter plane for the imperfect gluing 
bifurcation in the planar case. A one parameter family of 
(imperfect) systems, S, is indicated by a curve through the 
plane close to e = 0. The arrows indicate the direction in 
which orbits are created. 




FIG. 7: Bifurcation diagram (period against parameter) on 
the one-parameter path S of Figure ^ 



Jacobian matrix of the flow which are closest to the imag- 
inary axis. If, up to complex conjugation, these are Ai 
and A2 with Re Ai < < Re A2 then the saddle index, 6, 
defined by 



-Re Ai/Re A2 



(1) 



plays an important role. The two-parameter space near 
the intersection of Go and Gi in the planar case is shown 
in Figure ^ (Ai and A2 are real), where the symmetry 
is a point symmetry about the stationary point and the 
direction of time may be chosen so that ^ > 1. Each 
simple homoclinic bifurcation creates a periodic orbit in 
the direction indicated by the arrow on the bifurcation 
curve. The parameter plane is divided into six regions by 
the curves of bifurcations, and the periodic orbits (from 
the local theory) which exist in each region are indicated 
by their codes. The bifurcations observed on the one- 
parameter path S in Figure ^ are shown in Figure |^, 
which is the more conventional representation. 
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FIG. 8: Bifurcation diagram (period against parameter) of 
the modified global bifurcation as suggested by Figure |2| 



B. Relationship with the experiment 

The curves sketched in Figure are in reasonably good 
agreement with the experimental ones in Figure |^ except 
for the extra complications at larger parameter values 
described in section II D. Also the fact, mentioned at the 



end of section II C, that all three of the orbits labelled 
'0', '1' and '10' coexist for some values of ai. However, 
even these aspects can be incorporated into our picture 
of imperfect global bifurcations. For smaller values of /3i 
Shil'nikov wiggles are observed, suggesting that S < 1 
(and Ai is complex) in this parameter regime. In this 
case, as earlier, there may be symmetry-breaking and re- 
verse symmetry-breaking bifurcations of the symmetric 
orbit (in the perfectly symmetric system) |T7| . The bi- 
furcations observed in Figure and described in more 
detail in Figure ^ are not in the asymptotic region of ap- 
plicability of the homoclinic theory (large period, close to 
homoclinic bifurcation) and so we invoke an extra pair of 
assumptions on the underlying symmetric system for our 
model: that there is a symmetry-breaking and reverse- 
symmetry breaking bifurcation on the symmetric orbit 
and that S < 1. 

If (5 < 1 then the curves of homoclinic bifurcations are 
essentially as in Figure ^ but the direction of the bifurca- 
tions is reversed (more precisely, the diagram is reflected 
about the e— axis) and the orbits created are saddles 
(rather than stable, as would be the case ii 5 > 1). This 
now suggests the new interpretation of Figure ^ which is 
shown in Figure ^. The major new feature is that since 
the orbits are created in the opposite direction to the case 
with 5 > 1 in Figure ^ and are unstable, the points at 
which the orbits cannot be followed further (ai = 0.6041 
for the '10' orbit and ai = 0.6025 for the '0' and '1' orbits 
in Figure |^) are now assumed to be saddle-node bifurca- 
tions. There are a number of possible interpretations for 
the disconnected symmetry-breaking bifurcations, and 
one of these is shown in Figure although we make no 
claim that it is the most likely. Note that the new ar- 
rangement of the homoclinic bifurcations does provide a 
region of parameters where the orbits '0', '1', and '10' 
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coexist and are stable, as seen in the experiment. 

The important feature of the analysis above is that two 
assumptions about the underlying mathematical model 
are sufficient to explain the orbits observed in the ex- 
periment. It is worth emphasising that this can be done 
without constructing the model equations explicitly, sim- 
ply be suggesting that any model equation must have 
various dynamical features. 



C. Other cases 

In the literature, codimension two global bifurcations 
are generally described with Gq and Gi as the coordi- 
nate axes of the bifurcation analysis. In this case the 
symmetric system may be assumed to lie on the diagonal 
of the parameter space, with the asymmetry parameter 
perpendicular to the diagonal (just tilt the diagrams by 
45° to get an impression of the locus of bifurcations). It 
is, however, important to bear in mind that the curves 
Go and Gi appear to intersect with a very small angle 
of intersection in asymmetric perturbations of symmetric 
systems, whereas the standard analysis depicts the inter- 
section angle to be at right angles. Provided the intersec- 
tion is transversal the analysis is not affected, although it 
does mean that the true picture for the asymmetric per- 
turbation will be a very skewed version of the standard 
pictures. 

The basic feature common to all the relevant types of 
bifurcation we consider is that as the bifurcation param- 
eter /J, is varied, a (more or less complicated) sequence of 
bifurcation is observed with the net effect that a pair of 
periodic orbits (those we have labelled '0' and '1') is de- 
stroyed, and a single large periodic orbit is created. The 
precise details of the bifurcations depends on the system, 
but it is still possible to make a number of general state- 
ments. 



1. The one-sided case 

If the direction of time can be chosen so that A2 is real 
and 6 > 1 (cf. (1)) then the codimension one bifurcations 
on Go and Gi are one-sided and fairly general statements 
about the bifurcations involved in the range of validity of 
the rigorous argument: large period and parameters close 
to the intersection of Go and Gi are possible 22|. First, 
there are at most two periodic orbits, and second, any 
periodic orbit has a very particular description in terms 
of the symbols '0' and '1' introduced above. Technically, 
the sequences are rotation compatible sequences |Q , but 
in practice a simple consequence is that periodic orbits 
have codes of the form 



limit, p, of the number of Is in the sequence to the length 
of the sequence exists and is called the rotation number 
of the orbit. In one case (the so-called stable orientable 
Lorenz-like case, see [^), there is an infinite set of bi- 
furcations along a typical path and at any one parame- 
ter after crossing the first bifurcation curve, there is at 
most one periodic orbit. Moreover, the rotation number 
varies continuously along the bifurcation path, implying 
the existence of parameter values with non-periodic (but 
non-chaotic) attractors. 

If Ai is complex then the range of bifurcations is 
more complicated and depends on the precise path taken 
through the parameter space. Here there are regions 
of coexistence of certain periodic orbits - those whose 
rotation numbers — and — are Farey neighbours, i.e. 
\pi<l2 — 91P2I = 1 - but typical curves in parameter space 
do not intersect most of these regions. A more complete 
list of the possibilities can be found in |l| . 

All the bifurcations of the rigorous analysis involve 
one-sided global bifurcations, and there are no local bi- 
furcations on the branches of each periodic orbit. If these 
occur it is necessary to appeal to effects outside the rigor- 
ous region of validity of the mathematical results - this is 
made much easier by an understanding of the two-sided 
bifurcations. 



2. The two-sided case: Shil'nikov's wi 



Oini01"201»301»401^ 



(2) 



where for all i, G {m, m -\- 1} for some m > (or the 
same with the roles of and 1 exchanged). Moreover, the 



The symmetric bifurcation diagram of the Shil'nikov 
case (Ai complex, A2 real and (5 < 1) is given in 
The locus of the pair of orbits ('0' and '1') in parameter- 
period space oscillates as the period increases to infin- 
ity, with period-doubling and reverse period-doubling bi- 
furcations on every other branch. The symmetric orbit 
oscillates in a similar way, but with symmetry-breaking 
bifurcations on every other branch. Breaking the sym- 
metry of the system will have two effects ~ the global 
bifurcations which coincide in the symmetric system will 
be split apart and the symmetry-breaking bifurcations 
will typically become disconnected as described above. 
In the two-parameter diagram close to the intersection 
of Go and Gi, curves of more complicated bifurcations 
(Goi and Gio) oscillate rapidly and intersect each other 
(there are infinitely many other curves of homoclinic bi- 
furcations to complicate matters further). For a typical 
asymmetric path there will be a single intersection with 
Go and Gi , but potentially several intersections with Gio 
and Goi. The orbits created in the bifurcations involving 
Go and Gi will lie on the usual Shil'nikov wiggle in the 
parameter-period plane as observed experimentally (see 
Figure The symmetric orbit, '10', can also be fol- 
lowed experimentally (see Figure [l0| ); there are multiple 
intersections of the parameter path with Giq, i.e. extra 
bifurcations which create and destroy the orbits labelled 
'10*'. Between the conjectured intersection of the pa- 
rameter path with Go and Gio it is possible to observe a 
stable orbit with code '100'. Such an orbit can be created 
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FIG. 9: Coexisiting Shil'nikov wiggle at /3 = 0.5317. The 
branch noted with [ ] correspond to the '1' and with ( ) to the 
'0' orbit respectively. 
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FIG. 10; Shil'nikov wiggle and gluing process of the '10' orbit 
at /3 = 0.5317. The period of the '10' orbits and the '100' 
orbit is rescaled by two and three respectively. 
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FIG. 11: Time series and phase portraits of different dy- 
namical states involved in imperfect gluing bifurcation at 
(3 = 0.5317. (a) periodic orbit '0' on asymmetric branch at 
Qfi = 0.56119, (b) chaos at ai = 0.56125, (c) period-3 or- 
bit '100' at Qi = 0.56136, (d) chaos at qi = 0.56146, (e) 
symmetric periodic orbit '10' at qi = 0.56152 



from homoclinic orbits obtained from the gluing of the 
orbits '10' and '0'. These bifurcations are expected due 
to the intersection of Go and Gio in the two parameter 
analysis (cf. the S > 1 case in [^) which create extra 
curves of homoclinic orbits Goio and Gioo- 



IV. EXPERIMENT II 

The gluing process in the non-planar region of param- 
eter space involves complicated orbits which evolve on 
a Shil'nikov wiggle. A pair of these wiggles are shown 
in Figure ^ where the period is plotted as a function of 
ai at fixed /3i — 0.5317. Here the period of the orbit 
approaches infinity through a sequence of folds where al- 
ternate branches are unstable and indicated by dashed 
schematic lines in the figure. The stable solutions un- 
dergo forward and reverse period-doubling sequences on 
the first two folds whereas the highest period orbits only 



exist over a tiny range of the parameter. In a perfectly 
symmetric system these two wiggles would overlap com- 
pletely. The effect of the imperfection in the circuit is to 
displace the two curves from one another. 

A Shil'nikov wiggle has also been observed on the sym- 
metric orbit and the results are shown in Figure 
There we can see three levels of the wiggle with period 
doubling sequences. The '10' orbits in this case were 
asymmetric but we were unable to find the mirror image 
pairs of solutions in this case. We were, however, able to 
observe them at smaller values of Pi . The gluing process 
takes place on the third level with intervening sequences 
of chaos and a stable '100' or bit; a s expected from the dis- 
cussion at the end of section III C , Note we also observed 
the '10*' which is an i ntegra l part of the gluing process 
as discussed in section [II C above. A set of time-series 
and phase portraits are displayed in Figure |ll|. The '0' 
orbit on the disconnected branch glues to the '10' large 
scale orbit via two chaotic phases with an intermediate 
period-3 '100' sequence. 
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V. CONCLUSION 

Although symmetric equations are frequently used to 
model almost symmetric systems, we have shown that a 
more careful examination of experiments can reveal fea- 
tures which do not appear in the symmetric models. In 
particular, we have focussed here on global bifurcations 
which involve periodic states of the system, and we have 
shown how a number of complicated bifurcation diagrams 
observed in the experiments can be interpreted by ap- 
pealing to a theory of imperfect homoclinic bifurcations. 

A standard approach to the modelling of physical phe- 
nomena is to construct a mathematical model of the ex- 
periment, and use this to either predict or explain fea- 
tures of the experiment. This entails both the construc- 
tion of the model and the analysis of the model con- 



structed. It is noticeable that in the approach taken here 
we have appealed to properties of a model without having 
to either construct or analyse the model. We have sim- 
ply said that any mathematical model of the experiments 
must have certain features, and that these features lead 
to certain conclusions by the application of global bifur- 
cation theory. 



tion 



Bi furc ation diagrams consistent with those of sec- 
have now been observed in more physically in- 
teresting systems. Abshagen has found bifurcation 
diagrams with a striking similarity to Figure ^ in exper- 
imental data from fluid flow. We believe that the ap- 
proach taken here will find application in a broad variety 
of experiments in which symmetry, or rather, almost- 
symmetry, plays a role. 
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